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ABSTRACT 

We investigate the Sloan Digital Sky Survey (SDSS) photometry from Data Release 8 (DR8) 
in the search for systematic trends that still exist after the calibration effort of Padmanabhan et 
al. We consider both the aperture and point-spread function (PSF) magnitudes in DR8. Using 
the objects with repeat observations, we find that a large proportion of the aperture magnitudes 
suffer a ^0.2-2% systematic trend as a function of PSF full-width half-maximum (FWHM), 
the amplitude of which increases for fainter objects. Analysis of the PSF magnitudes reveals 
more complicated systematic trends of similar amplitude as a function of PSF FWHM and 
object brightness. We suspect that sky over-subtraction is the cause of the largest amplitude 
trends as a function of PSF FWHM. We also detect systematic trends as a function of sub- 
pixel coordinates for the PSF magnitudes with peak-to-peak amplitudes of ^1.6 mmag and 
~4-7 mmag for the over- and under-sampled images, respectively. We note that the systematic 
trends are similar in amplitude to the reported ~1% and ^2% precision of the SDSS photom- 
etry in the griz and u wavebands, respectively, and therefore their correction has the potential 
to substantially improve the SDSS photometric precision. We provide an IDL program specif- 
ically for this purpose. Finally, we note that the SDSS aperture and PSF magnitude scales are 
related by a non-linear transformation that departs from linearity by ~l-4%, which, without 
correction, invalidates the application of a photometric calibration model derived from the 
aperture magnitudes to the PSF magnitudes, as has been done for SDSS DR8. 

Key words: instrumentation: detectors - techniques: photometric - surveys. 



1 INTRODUCTION 

The SDSS dYork et ail |2000|) provides photometry (including 
widely-used aperture and PSF magnitudes) for objects down to 
r ~22.5 mag and covering 14555 deg 2 jAihara et al.l l201ll - 
from now on AIH1 1). The photometric calibration is uniform over 
the survey area at the ~1% and ~2% levels in t he griz and u 
wavebands, respectively ([Padmanabhan et al. 2008 - from now on 
PAD08). This vast data set has provided insight into many areas 
of astronomical research, including the large-scale structure of the 
Universe, properties of galaxies, Galactic structure, stellar popula- 
tions, just to name a few. 

For many of the scientific applications of SDSS photometric 
data, an accurate absolute photometric calibration (i.e. knowledge 
of fluxes in physical units such as J raT 2 s _1 ) is not as impor- 
tant as a precise relative calibration over the survey. Furthermore, 
it is substantially more challenging to achieve an accurate abso- 
lute calibration due to the inherent difficulties in minimising the 
systematic uncertainties in the mapping between the spectral en- 
ergy distribution (SED) of an appropriate (and usually very bright) 
fundamental spectrophotometric standard and the natural photo- 
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metric system of the survey. For these reasons, much recent ef- 
fort has been invested in determining a precise relative calibra- 
tion of SDSS i maging data, and ot her similar surveys have fol- 
lowed suit (e.g. lRegnault et al . 2009). Future ground-based all-sky 
surveys that require 1% or better relati ve photo metric calibrations 
(e.g. PanSTARRS - iKaiser et alj|2002l . LSST - llvezic et al]|2008t 
etc.) are investing considerable resources to perform exquisitely 
precise measurements of the telescope-instrument-detector system 
throughput as a function of wavelength for each imager pixel 
( Stubb s & Toru^2006tlSmbbs" et al.l2010l) and to use spectroscopic 
observations of calibrator stars coupled with detailed atmospheric 
modelling to enable the precise correction of atmosp heric extinc- 
tion effects to the leve l of a few millimagnitudes ( Stub bs et al.l 
l2007l : lBurkeetal.ll2010l) . 

The purpose of a relative photometric calibration is to model 
the systematic differences in the measured magnitudes of survey 
objects as they are observed over time and over the spatial extent of 
the survey. The sources of these differences are multitude although 
they are usually dominated by atmospheric extinction variations, 
differences in detector responses, and errors introduced during the 
photometric analysis of the images. Ideally, correction of the full 
survey data using the fitted calibration model should result in ob- 
ject photometry that reaches some target homogeneity such as the 
desirable 1% level. The key to being able to fit a relative photomet- 
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ric calibration across survey data is the existence of repeat observa- 
tions of a sufficient number of objects to link all the observational 
data via the adopted calibration model and to provide strong enough 
constraints on the model parameters. The early development of the 
photometric calibrat ion methodology may be fo und in the works of 
ikarris et all dl98ll). iReed & Fitzgerakj Jl982l) [Manfroid & Heckl 
d 1 9831) . iHonevcuttU 1 992h . and lManfroid 1 1 995h . 

PAD08 present the impressive development of a photometric 
calibration model for SDSS data on a scale many times larger than 
has been used before. The model is fit to the repeat observations 
from the overlaps between the scan columns and at the ends of 
the scans along great circles (see Section 13. U , and from the re- 
peat scanning of certain scan re gions. For SDSS Data Release 6 
dAdelman-McCarfhv et ai]|2008l) . their model has ~2000 calibra- 
tion parameters of interest and ~10 7 nuisance parameters (the star 
mean magnitudes) to be solved for. PAD08 marginalise over the 
unknown star mean magnitudes and iteratively solve the normal 
equations of the smaller least squares problem for the calibration 
parameters. The resulting relative calibration is good to ~1% in 
the griz wavebands (~2% in the u waveband). 

In PAD08, the authors looked at the calibration residuals as a 
function of star brightness (or mean magnitude), as a function of 
detector column to assess the accuracy of the flat-field calibrations, 
as a function of time, and as a function of celestial coordinates, 
the latter being further decomposed into spatial error modes. They 
found a clear trend as a function of time, which manifests as "coher- 
ent errors at the few millimagnitude level", and which most likely 
correlates with time-varying image quality parameters such as see- 
ing, sky brightness, etc. An example of such a time-dependent trend 
may be discerned as a ~ 10 mmag drift in the magnitude residuals 
in the upper panel of Fig. 9 in PAD08. This implies that the PAD08 
photometric model may need to be updated to include appropriate 
terms to fit these trends. We also highlight the fact that the pho- 
tometric calibration parameters in PAD08 are derived from fitting 
the SDSS aperture magnitudes and the calibration is then applied 
to calibrate various object magnitude measurements in the survey 
catalogues such as PSF magnitudes. 

The purpose of this paper is to extend the investigation into 
the systematic trends in the SDSS photometric data. We stress that 
an exhaustive study of this topic is virtually impossible and not 
within the scope of our work. We start in Section |2]by describing 
our modelling method. In Section[3] we describe the procedure we 
use to select appropriate photometric data from the SDSS data base. 
We then investigate the systematic trends in the SDSS aperture and 
PSF magnitude measurements, and we report on the trends we find 
as a function of PSF FWHM, object brightness, and subpixel co- 
ordinates (Section|4j. In Section|5] we derive the relation between 
the aperture and PSF magnitude scales, which we find suffers from 
non-linearities. We summarise and discuss our results in Section[6] 



Instead, we introduce a binning for an object/image property 
X that covers the full range of data values that we are modelling. 
For the fcth bin in X, we introduce an unknown magnitude off- 
set Zk to be determined, the purpose of which is to model the mean 
difference of the photometric measurements within the correspond- 
ing bin in X from the rest of the photometric measurement^ Our 
photometric model may then be written as: 



rrii — S~] Sjv M p + 2~] dkp Z p = Mj + Z k 



(1) 



p=i 



p=i 



where rrU is the model magnitude for the ith magnitude measure- 
ment rrii, Mj is the unknown true instrumental magnitude of the 
jth object, A^bj is the number of objects in the photometric data 
sample, N z is the number of magnitude offsets Z^ that we are at- 
tempting to fit, and Sij is the Kronecker delta-function: 



1 if i — j 
if Mj 



(2) 



Note that the ith photometric measurement in our photometric data 
sample belongs to the j(i)fh object and the fc(i)fh bin in X, where 
the adopted notation for j and k reflects the fact that both of the 
indices j and k are functions of the index i. However, in the rest 
of this paper, we devolve to using the notation j and k for 
and k(i), respectively, in order to avoid confusion in our subscript 
notation. 

The photometric model defined in Equation [TJ is degenerate 
since increasing each of the instrumental magnitudes Mj by an ar- 
bitrary amount c at the same time as decreasing each of the mag- 
nitude offsets Zk also by c will have no effect on the model values 
rrii. This degeneracy may be removed by fixing the instrumental 
magnitude of a single object to an arbitrary but convenient value. 
This is equivalent to fixing the absolute photometric calibration of 
the observations. 



2.2 Constructing The Normal Equations 

Equation [T]is linear and th erefore we may use general linear least- 
squares flPress et al.ll20071) to find the solution for the instrumental 
magnitudes Mj and the magnitude offsets Z%. The chi-squared as- 
sociated with EquationQ]is: 



E 



J2p=l $kp Zp 



01 



(3) 



where iVdata is the number of photometric measurements, and ai is 
the uncertainty on the ith photometric measurement. 

The normal equations of the least-squares problem are found 
by differentiating the \ 2 in Equation[3]with respect to each param- 
eter in the model and they may be written in matrix form as: 



2 MODELLING THE PHOTOMETRIC DATA 
2.1 The Photometric Model 

Our approach to identifying the systematic trends in the SDSS pho- 
tometric data consists of investigating how the photometry depends 
on each object or image property which we suspect could have an 
impact on the photometry. Our aim is to identify any significant 
trends that were not modelled during the PAD08 calibration of the 
aperture magnitudes. However, since we do not know a priori the 
form of any potential trend, we cannot assume a smooth functional 
form at this stage. 



A 
B T 



(4) 



where A, B, and D are matrices of sizes ATobj x JV bj, iVobj X N z , and 
N z x N z elements, respectively, and where vi and V2 are vectors of 
lengths ATobj and AT Z elements, respectively. The vector xi contains 
the ATobj unknown instrumental magnitudes Mj , and the vector X2 



1 Note that a positive magnitude offset indicates that the photometric 
measurements within the corresponding bin in X are fainter on average 
than the rest of the photometric measurements, and vice versa. 
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contains the N z unknown magnitude offsets Zk ■ The individual el- 
ements of A, B, D, vi, and v 2 are given by: 





= S 1P 5 jq 1 °i 
i=l 


(5) 


B pq 


i=l 


(6) 


D pq 


= Skp Sk q / erf 

! = 1 


(7) 


V 1,P 


= ^ Sjp TUi 1 <j\ 
! = 1 


(8) 


V2,p 


= ^2 Skp rrii 1 al 

i=l 


(9) 



Matrix A (Equation[5} is diagonal since the product Sj P Sj q is 
only non-zero when j = p — q, which makes sense because a pho- 
tometric measurement m* may only correspond to a single object. 
Similarly, matrix D (Equation [7} is diagonal, which is due to the 
fact that a photometric measurement rrii may only correspond to a 
single bin in X. Matrix B (Equation [6J happens to be sparse when 
there are generally less photometric measurements per object than 
bins in X. Finally, we observe that the elements of the vectors vi 
and V2 (Equations[8]&[9} are simply the inverse-variance weighted 
sums of the photometric measurements rrii for each object and for 
each bin in X, respectively. 

At this stage, it is worth noting that the construction of the 
least squares matrix and the vector on the right-hand side is a highly 
parallel computational problem due to the natural partitioning of 
the photometric measurements on a per-object and per-bin-in-X 
basis, which is a very convenient fact because this is where a large 
fraction of the computational operations occur in finding the solu- 
tion for xi and Xa. Furthermore, the facts that A and D are diag- 
onal, and B is sparse, may be used to implement computer code 
that makes efficient use of available memory resources, especially 
when it may be unfeasible to store the full versions of the matrices 
A, B, and D in computer memory. 



2.3 Solving The Normal Equations 

The fact that A is diagonal in the normal equations for photomet- 
ric models including the true instrumen tal magnitudes ( as in Equa- 
tion[T} has been spott ed previous l y by [ Regnault et al] d2009l) . and 
it is clearly the case in Honeycutt ( 1992) even though no comment 
was made. In Appendix A of iRegnault et alj d2009l) . they use this 
property to develop a tractable way of solving the normal equa- 
tions for the parameter vectors xi and X2, and for calculating the 
marginalised covariance matrix C for the parameters in Xz. We 
briefly repeat this method here using the notation in this paper. 

Eliminating the parameter vector xi from the normal equa- 
tions in Equation|4]yields: 



D — B T A *B 



x 2 = v 2 



Vi 



(10) 



where the elements of A 1 are easily computed as: 



A" 1 - 

f-pq 



I /A pq ifp = q 
ifp ^ q 



Cholesky factorisation of the symmetric and positive-definit^fl 
matrix (D — B T A _1 B), followed by forward and back sub- 
stitution is the most ef ficient and numerically stable method 
dGolub & Van Loan! 1996b for obtaining the solution for the magni- 
tude offsets Zk ■ The solution for the instrumental magnitudes Mj , 
if required, can be obtained by substituting the solution for x 2 into 
Equation [4] 



Xi 



: A Vl 



A Bx 2 



(12) 



Finally, the marginalised covariance matrix C for the magnitude 
offsets Zk may be obtained by marginalising over the true instru- 
mental magnitudes: 



D 



B T A 



(13) 



2.4 Implementation And Processing Time 



We have implemented an IDL program 

f it_photometric_calibration . pro as part 



called 
of the 

DanIDL0 library of routines. This program allows the user to fit 
a fully configurable and highly flexible photometric calibration 
model (including s tatic and rotating illumination corrections - 
Moehler et alJl201Cl) to a set of magnitude measurements, and it 
employs the methodology developed in Sections 12.21 & 12.31 The 
specific photometric model defined in Equation[TJuses only one of 
a slew of available terms that may be configured in the photometric 
model in f it_photometric_calibration . pro. 

The program that we have developed works solely by storing 
all the necessary data and arrays in computer memory, which lim- 
its the size of the fitting problem that can be tackled to the amount 
of memory that is available. If we improved our program to capi- 
talise on the highly parallel nature of the computational problem by 
reading and writing data from the computer hard disk as necessary, 
then the limit on the size of the fitting problem that can be solved 
can be greatly increased to match with the amount of disk space 
available since computer hard disks currently obtain much larger 
storage sizes than computer memory. By doing this, we would also 
remove any constraint imposed by the amount of computer mem- 
ory available, because the program would only use a limited and 
well-defined amount of memory at any one time. However, for the 
analysis of the data presented in this paper, our IDL program is not 
limited by the amount of computer memory that is available. 

We have briefly tested the performance of our program and 
its scalability to larger data sets. Using a single Intel Xeon 2.0Ghz 
CPU on a 64-bit machine with 132Gb of RAM, we find that for 
N z — 75 and iVdata = 3iVobj, our program uses ~1.7, 4.4, and 
8.9 Gb of RAM for N obl = 0.6 x 10 6 , 1.5 x 10 6 , and 3 x 10 6 , respec- 
tively, while taking ~1.6, 3.8, and 8.5 minutes to run, respectively. 
The processing limit for our program on the same machine seems 
to lie at approximately jV ob j = 10 7 , iVdata = 10 s , and N z = 100, 
which requires -~48 Gb of RAM and takes ~52 minutes to run. 



3 SDSS PHOTOMETRIC DATA 
3.1 The SDSS Imaging Camera 

The SDSS imaging camera, mounted on a dedicated 2.5m telescope 
at the Apache Point Observatory (New Mexico), consists of an 



(11) 2 We do not provide the proof of this statement in this paper. 

3 http://www.danidl.co.uk 
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array of 30 SITe/Tektronix charge-coupled devices (CCDs), each 
of size 2048x2048 pix and with a pixel scale of 0.396" pix -1 , 
arranged in six columns of five chips each wi th a space of ap- 
proximately one chip width between columns dGunn et all 1 19981 : 
iGunn et al.l2 006). Each row of six chips is positioned behind a dif- 
ferent filter so that SDSS imagin g data are produced in five wave- 
bands , namely u, g, r, i, and z dFukugita et al.lfl996l ISmith at al.l 
120021) . The camera operates in time-delay-and-integrate (TDI) 
reado ut mode scanning along great circles in the sky at the sidereal 
rate dZaritskv. Shectman & B redthauer 1996). The camera scan di- 
rection is parallel to the six columns of CCDs, and the temporal 
observation order of the wavebands is riuzg. The chip arrange- 
ment is such that two scans may be used to cover a filled stripe 
2? 54 wide, with ~l'(~8%) overlap between chip columns in the 
two scans. The telescope and camera optics, along with the CCD 
layout, are designed such that, given perfect tracking along great 
circles, the optical distortion over the field-of-view (FOV) results 
in star tracks across each CCD that deviate from being parallel to 
the CCD columns by at most ~0.06"(or ~0. 15 pix) over the length 
of a detector in the worst case. 



3.2 Data Selection From SDSS DR8 

We obtained our photometric data sample from SDSS DR8 
(AIH11), which includes all imaging data up to the retire- 
ment of the SDSS camera, reprocessed with an updated ver- 
sion of the SDSS photometric pipeline, and calibrated using the 
PAD08 modelling scheme. We employed the CasJobs interfac^B 
to query the DR8 data base. We extracted photometric obser- 
vations of all starfl (TYPE=6) th at have no "child objects" 
(NCHILD=0; IStoughton etai]|2002h and calibrated (CALIB STA- 
TUS .FILTER AND 1) PSF magnitudes brighter than 19 mag 
(PSFMAG_FILTER<19) with positive uncertainties smaller than 
1 mag (0<PSFMAGERR_FILTER<1). The magnitude threshold 
was chosen to limit our data sample to those objects with photo- 
metric uncertainties of better than ~2-3%. The photometric mea- 
surements are supplied with corresponding quality flags which aid 
in the selection of "good quality" measurements. We further fil- 
tered our dat a by applying the set o f quality constraints described 
in Table 2 of Brami ch et alj d2008l) . Finally, the resolution of the 
individual objects to which each photometric observation belongs 
is available in DR8 and the object identifications are stored in the 
THINGID entry in the DR8 data base. Where this identification 
process fails, a THINGID value of — 1 is stored, and we therefore 
dropped such measurements from our data sample. 

The aperture photometry measurements in the SDSS DR8 data 
base are stored as fluxes in units of "nanomaggies", defined such 
that an object of brightness 1 nanomaggie has a magnitude of 22.5 
in any waveband. In our data sample, we included the aperture pho- 
tometry measurements for the commonly used "aperture-7" aper- 
ture of radius 7.43", and we converted these flux measurements and 
their uncertainti es into standard SDSS " a sinh m agnitudes" using 
the equations in iLupton. Gunn & Szalavl (l999) and the standard 
SDSS values for the "softening" parameter 6 of 0.14, 0.09, 0.12, 
0.18, and 0.74 nanomaggies in ugriz, respectively. We note that 
our data are selected on the basis of PSF magnitudes, and some- 
times the SDSS pipeline failed to measure a corresponding aper- 

4 http://skyservice.pha.jhu.edu/casjobs/ 

5 In reality, we are selecting PSF-like point-source objects which may in- 
clude quasars, asteroids, etc. 



ture flux. We flagged these cases for later exclusion when we wish 
to perform calculations using the aperture magnitudes. 

We also performed a neighbourhood search within our data for 
each waveband independently, recording for each object the num- 
ber of objects within a radius of 15", and storing this information 
with each photometric observation in our data set. 

Our final data set consists of 30 data tables, one for each detec- 
tor, listing for each object observation the object identification, the 
right ascension and declination, the aperture and PSF magnitude 
measurements and their uncertainties, the PSF FWHM, the detec- 
tor x and y coordinates, and the number of objects from our data 
set within 15". In Table[T] we report the number of objects and the 
number of associated photometric observations in our data set for 
each detector. These quantities differ between the PSF and aperture 
magnitude measurements due to the failure of the SDSS pipeline to 
measure an aperture flux in some cases. 

In Fig. I l(a)[ we pl°t a histogram of the PSF FWHM (arcsec) 
over the object observations for each filter. The vertical dashed lines 
represent the median PSF FWHMs of ~ 1.46", 1.36", 1.24", 1.18", 
and 1.20" for the ugriz wavebands, respectively. We note that the 
actual PSF FWHM distributions for each detector among detectors 
with the same filter differ at the ~5-20% level due to PSF varia- 
tions across the camera FOV, with the detectors in the outer camera 
columns having a larger PSF FWHM than those in the inner cam- 
era columns. This may be seen in columns 7, 8, and 9 of Table [T] 
where we list the 5th, 50th (median), and 95th percentile values of 
the cumulative distribution of PSF FWHM values for each detector, 
respectively. 

In Fig. |l(b)| we plot a histogram of the number of observations 
of each object over the set of objects for each filter. Clearly most 
objects have only a single observation in SDSS DR8. Those objects 
that do have repeat observations generally have only two or three 
repeat observations, highlighting the fact that in SDSS DR8, it is the 
overlaps between scan columns and at the ends of the scans along 
great circles that dominate the generation of repeat observations. 

3.3 Repeat Observations 

In the absence of external information about the relative bright- 
nesses of objects, it is only possible to analyse the systematic trends 
present in a set of photometric measurements by using those objects 
for which repeat observations have been performed. An object with 
a single photometric observation simply provides one data point 
at the expense of generating one unknown (the object brightness), 
bringing no new information to the problem. Hence, for our analy- 
sis of the systematic trends present in the SDSS aperture and PSF 
magnitudes, we must select the subset of repeat observations from 
our full DR8 data set. 

We work with the data for each detector independently. For 
the aperture magnitude measurements, we select the observations 
for which the SDSS pipeline successfully measured an aperture 
flux, and then further restrict our selection to those observations 
for which an object has at least two observations. We again pro- 
duce 30 data tables, one for each detector, and in Table|2] we report 
the number of objects and the number of associated photometric 
observations in this data set for each detector. 

For the PSF magnitudes, again we select the observations for 
which an object has at least two observations, and we produce an- 
other 30 data tables, one for each detector. In Table [3] we report 
the number of objects and the number of associated photometric 
observations in this data set for each detector. 

We note that the PSF FWHM distribution of the repeat obser- 
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Figure 1. (a) Histograms of PSF FWHM (arcsec) over the object observations for each filter, using data for the six camera columns combined. The vertical 
dashed lines represent the median PSF FWHMs for each filter, and the vertical dotted black line represents the PSF FWHM corresponding to critical sampling 
(0.93") for the SDSS imaging detectors, (b) Histograms of the number of observations of each object over the set of objects for each filter, using data for the 
six camera columns combined. Both: The data used in these plots are the observations we selected from SDSS DR8 in Section |3~2l 



Table 1. The properties of the selected observations from SDSS DR8 (see Section |3~2l for each of the 30 SDSS imaging detectors. The data sets for each 
detector are organised by filter and camera column. The 5th, 50th (median), and 95th percentile values of the cumulative distribution of PSF FWHM values 
(arcsec) are listed in columns 7, 8, and 9, respectively. 



Filter Camera PSF Magnitudes Aperture Magnitudes PSF FWHM (arcsec) 

Column No. Of Objects No. Of Observations No. Of Objects No. Of Observations 5th Percentile Median 95th Percentile 



u 


1 


2552954 


3680087 


2540345 


3661050 


1.21 


1.56 


2.19 


u 


2 


2717648 


3970893 


2698799 


3941161 


1.17 


1.50 


2.15 


u 


3 


2694623 


3944556 


2678855 


3920169 


1.05 


1.40 


2.07 


u 


4 


2743806 


4020514 


2726917 


3994001 


1.04 


1.39 


2.05 


u 


5 


2719071 


3964543 


2702696 


3939132 


1.05 


1.40 


2.05 


u 


6 


2786132 


4092273 


2768708 


4064715 


1.13 


1.47 


2.10 


9 


1 


5511948 


7946945 


5432614 


7835474 


1.07 


1.39 


2.00 


9 


2 


5515143 


7978389 


5429156 


7856395 


1.06 


1.36 


1.97 


9 


3 


5457676 


7889987 


5383166 


7783855 


1.03 


1.34 


1.94 


9 


4 


5434436 


7859519 


5360404 


7753636 


1.02 


1.33 


1.93 


9 


5 


5632190 


8169352 


5552972 


8056549 


1.03 


1.35 


1.95 


9 


6 


5608587 


8139232 


5535674 


8036950 


1.06 


1.38 


1.97 


r 


1 


8142785 


11698216 


7492919 


10718162 


0.94 


1.27 


1.83 


r 


2 


8061099 


11595390 


7377343 


10558970 


0.90 


1.23 


1.79 


r 


3 


8067250 


11635364 


7384481 


10558611 


0.87 


1.18 


1.76 


r 


4 


8103852 


11695080 


7434713 


10668356 


0.88 


1.19 


1.75 


r 


5 


7761595 


11210722 


7092689 


10191515 


0.93 


1.24 


1.80 


r 


6 


8409320 


12184076 


7770240 


11218480 


1.01 


1.34 


1.87 




1 


10424397 


14944150 


9689614 


13844569 


0.88 


1.21 


1.75 




2 


10737931 


15545164 


9959987 


14360505 


0.82 


1.13 


1.68 




3 


10198773 


14609342 


9484867 


13532858 


0.80 


1.10 


1.65 




4 


10236481 


14778349 


9519267 


13683502 


0.81 


1.12 


1.66 




5 


1 1060536 


16060059 


10274455 


14865032 


0.86 


1.16 


1.70 


i 


6 


10712588 


15516036 


9991929 


14435648 


0.99 


1.33 


1.84 


z 


1 


13596773 


19575551 


11964994 


17056492 


0.96 


1.27 


1.77 


z 


2 


13679274 
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17082888 


0.88 


1.17 


1.68 


z 


3 


13745026 


19896381 


11997206 


17169875 


0.83 


1.11 


1.62 


z 


4 


13922390 


20149555 


12046553 


17227888 


0.85 


1.13 


1.63 


z 


5 


13593815 


19575105 


11946407 


17040156 


0.88 


1.18 


1.68 


z 


6 


13359060 


19139276 


11765172 


16715773 


0.97 


1.32 


1.80 
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vations is very similar to the PSF FWHM distribution of all selected 
DR8 observations displayed in Fig. |l(a)| which means that by us- 
ing the subset of repeat observations to investigate the systematic 
trends in the SDSS photometric data, we are not introducing a bias 
into our results. 



4 SYSTEMATIC PHOTOMETRIC TRENDS 
4.1 Aperture Magnitudes 

We may now investigate how the DR8 aperture photometry corre- 
lates with various object/image properties to see if there are any sig- 
nificant systematic trends that were not modelled during the PAD08 
calibration of the aperture magnitudes. Even though the aperture 
magnitudes themselves are rarely used for high precision photom- 
etry, this investigation is crucial because the fit of the PAD08 cal- 
ibration model to the aperture magnitudes is used to calibrate var- 
ious object magnitudes, including the PSF magnitudes, which are 
widely used for high precision photometry. We note that we cannot 
investigate the existence of systematic trends as a function of object 
brightness since we are using the photometry of objects of suppos- 
edly constant brightness, and such observations do not contain any 
information on the magnitude scale itself. In other words, if there 
are any non-linearities in the aperture magnitude scale, we cannot 
detect them using the repeat observations. 

We use the data set corresponding to Table [2] and we anal- 
yse the photometric data for each detector independently. For each 
quantity X that we wish to investigate, we fit the photometric 
model described in Section |2~T| to the aperture magnitude measure- 
ments using an iterative procedure where we reject photometric ob- 
servations that lie more than 3a away from the fitted model. We 
then drop the object observations for which an object now has only 
a single observation. This data rejection step is necessary because 
of the presence of variable sources and outlier photometric mea- 
surements in our data set, and typically ~3-5% of the data are re- 
jected at this stage. We then repeat the fit of the photometric model 
to the cleaned aperture magnitude measurements. 

Applying this iterative fitting procedure to the quantity X 
representing the detector row with bins of size 10 pix, we find 
that there are no significant trends in the fitted magnitude offsets 
(Zk < 0.2 mmag). Also, for the quantity X representing the detec- 
tor column with bins of size 10 pix, we find that there are no signif- 
icant trends in the magnitude offsets (Zk < 0.5 mmag). However, 
we do find systematic trends in the magnitude offsets when we per- 
form the fit for the quantity X representing the PSF FWHM. 



4.1.1 Trends As A Function Of PSF FWHM 

For each detector, we partition our data into bins of width 0. 1" in 
PSF FWHM, and we further subdivide the data in each bin into 
those observations which have magnitude measurements brighter 
than 16 mag, in the range 16-18 mag, and fainter than 18 mag. 
We construct our photometric model as in Section |2~T1 adoptin g an 
unknown magnitude offset for each bin, and we follow the iterated 
fitting procedure described in Section |4~T1 

In Fig. [2] for camera column 3, we plot the fitted magnitude 
offsets for the magnitude measurements brighter than 16 mag as 
a function of PSF FWHM (black points). We also plot the fitted 
magnitude offsets for the magnitude measurements in the range 
16-18 mag (red points), and fainter than 18 mag (green points). 



In each plot, the error bars represent the uncertainties in the magni- 
tude offsets, and they are generally much smaller than the plot sym- 
bols (~0.02-0. 1 mmag). We refrain from plotting magnitude offsets 
with uncertainties of greater than 5 mmag, and we mark those mag- 
nitude offsets that fall outside of the plot range with an asterisk of 
the relevant colour. The vertical blue line in each plot represents the 
PSF FWHM corresponding to critical sampling (0.93")- Objects in 
the sky scan from left to right in our plot layout. The corresponding 
plots for the remaining camera columns are very similar, and for 
brevity, we do not reproduce them here. 

We find in general that as the PSF FWHM increases, the mag- 
nitude offsets increase, indicating that the larger the PSF FWHM, 
the fainter an object is measured by the SDSS pipeline aper- 
ture photometry routines with the 7.43" fixed aperture. We also 
find that this trend is much stronger for fainter objects, with the 
strongest trend exisiting in the z waveband. The amplitude of 
this trend over the PSF FWHM range is ~7-15 mmag in each 
waveband for the magnitude measurements brighter than 16 mag, 
~ 12-25 mmag in each waveband for the magnitude measurements 
in the range 16-18 mag, ~30-60 mmag in the ugri wavebands 
for the magnitude measurements fainter than 18 mag, and ~100- 
170 mmag in the z waveband for the magnitude measurements 
fainter than 18 mag. We have also tested that the form and ampli- 
tude of this trend is independent of right ascension, declination, de- 
tector coordinates, subpixel coordinates, and the number of neigh- 
bouring objects within 15" by making similar plots to those dis- 
played in Fig. [2] for various mutually-exclusive partitions of our 
photometric data. 

The detected systematic trend in the magnitude offsets as 
a function of PSF FWHM is consistent with a problem of sky- 
background over-subtraction in the images from which the aperture 
fluxes are calculated. Sky over- subtraction by a constant amount 
per pixel in an image area containing objects of a range of bright- 
nesses implies that each of these objects will have an aperture 
flux measurement that is under-estimated by a constant amount 
A/, and therefore the ratio A///, where / is the true object 
aperture flux, will be larger for the fainter objects. The result is 
that the fainter objects will have a larger magnitude over-estimate 
Am ~ — 2.51og 10 (l — A///) than the brighter objects. As we 
have stressed in Section |4~TI we cannot detect such a distortion of 
the aperture magnitude scale itself using the repeat observations 
of objects of constant brightness. However, we can detect this ef- 
fect through the amplification of any PSF FWHM dependent trend 
in the magnitude offsets for fainter aperture magnitude measure- 
ments. In our case, the amplification of the trend manifests itself as 
a stronger gradient in the magnitude offsets as a function of PSF 
FWHM for the fainter magnitude measurements. 

The fact that objects are measured as systematically fainter 
as the PSF FWHM increases can also easily explained by a sky- 
background over- subtraction problem. An image with a larger PSF 
FWHM has the flux from objects spread out over more image pixels 
than an image with a smaller PSF FWHM, leading to a larger frac- 
tion of supposedly pristine sky-background pixels that are contam- 
inated by a small amount of flux from the PSF wings of the image 
objects. Inclusion of such contaminated pixels in the sky estimation 
is hard to avoid due to the very low level of object flux contamina- 
tion over a large area, resulting in a systematic over-estimate (and 
therefore over- subtraction) of the sky background level which is 
manifestedly worse for images with a larger PSF FWHM. 

AIH 1 1 discuss the fact that the sky estimation procedure of the 
SDSS pipeline has a tendency to over-subtract the outer regions of 
large galaxies on the sky as found by several authors (see references 
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Table 2. The properties of the selected repeat observations with aperture magnitude measurements from SDSS DR8 (see Section |3~3t for each of the 30 SDSS 
imaging detectors. The data sets for each detector are organised by filter and camera column. In columns 5-7, we list the 4<r upper limits L U nder, ^normal, and 
L wc ll that we determine on the existence of systematic photometric trends as a function of subpixel coordinates for the regimes of under-sampled, normally- 
sampled, and well-sampled PSFs, respectively. Upper limits are not reported for the r waveband because we actually detect a systematic trend in this case (see 
Section |4.1.2) . In column 8, we report the amount by which the \ 2 for the repeat observations decreases after calibrating the aperture magnitudes using the 
fitted magnitude offsets from Section l4. 1.11 In column 9, we list the A\ 2 values in units of a. 
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Figure 2. Systematic trends in the fitted magnitude offsets for the aperture magnitude measurements from camera column 3 as a function of PSF FWHM, 
object brightness, and subpixel coordinates (see text in Sections |4. 1 . 1 1 & |4. 1 .21 for details). The data used to generate these plots are the selected repeat 
observations with aperture magnitude measurements from SDSS DR8 (see Section l3~3t . The black, red, and green points correspond to aperture magnitude 
measurements brighter than 16 mag, in the range 16-18 mag, and fainter than 18 mag, respectively. Asterisks represent magnitude offsets that fall outside of 
the plot range. The three images above each plot represent the IPSMs for the regimes of under-sampled, normally-sampled, and well-sampled PSFs (from left 
to right). The colour-scale is linear and runs between =p4cr and =pl2cr for the ugiz and r wavebands, respectively. The full version of this figure including all 
of the camera columns is available from the authors on request. 
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therein). We suspect that the systematic trend that we have detected 
here is a symptom of the same problem for unresolved objects. 

We have also considered the possibility that the detected sys- 
tematic trend in the magnitude offsets as a function of PSF FWHM 
is due to increased flux losses outside of the fixed aperture when 
the PSF FWHM is larger. However, if this was the case, then the 
trend should be magnitude independent, since the fraction of ob- 
ject flux lost outside of the fixed aperture is independent of the ob- 
ject brightness. This is contrary to what we observe, and therefore 
we may be confident that this effect is not the principal cause of 
the detected systematic trend, although it may still contribute to the 
trend. Correction of the sky-background over-subtraction problem 
at the image processing stage will help to reveal whether aperture 
flux losses affect the aperture flux measurements as a function of 
PSF FWHM for the 7.43" fixed aperture. 

We also mention one other possible explanation for the de- 
tected systematic trend. Aperture flux losses for isolated objects are 
a function of PSF FWHM only. However, for objects that have one 
or more neighbouring objects whose fluxes contaminate the fixed 
aperture, the amount of flux contamination, and therefore the mea- 
sured aperture flux, will have a complicated dependence on PSF 
FWHM (the exact form of which is determined by the spatial dis- 
tribution of the neighbouring objects relative to the aperture bound- 
ary), and the aperture magnitude measurements for fainter objects 
of interest will be affected by a larger amount by the flux contam- 
ination. The fact that we have already found earlier in this section 
that the form and amplitude of the detected systematic trend in the 
fitted magnitude offsets as a function of PSF FWHM is indepen- 
dent of the number of neighbouring objects within 15" allows us 
to rule out this explanation as a primary cause. 

4.1.2 Trends As A Function Of Subpixel Coordinates 

Now we investigate the systematic trends in the aperture magni- 
tudes as a function of subpixel coordinates. For each detector, we 
partition our data into a uniform grid of 5 x 5 bins in subpixel coor- 
dinates covering the area of a single detector pixel using the object 
centroidfl and we further subdivide the data in each bin into those 
observations for which the corresponding PSF FWHM is less than 
the critical sampling of 0.93", greater than the critical sampling but 
less than 1.5", and greater than 1.5". We construct our photometric 
model as in Section l2~T"l adopting an unknown magnitude offset for 
each bin, and we follow the iterated fitting procedure described in 
Section|4~T1 

In Fig. [2] above each plot of magnitude offsets versus PSF 
FWHM, there are three square "images" with 5x5 elements. Each 
image represents the area of a single detector pixel, orientated such 
that objects in the sky, which scan along the detector columns, scan 
from left to right. We display the fitted magnitude offsets for the ob- 
servations with a PSF FWHM that is less than the critical sampling, 
greater than the critical sampling but less than 1.5", and greater 
than 1.5", in the left-hand, middle, and right-hand images, respec- 
tively. The image values are displayed in units of a, where o is 

6 We calculate the subpixel coordinates of the object centroid as the frac- 
tional part of the pixel coordinates of the object centroid, which are stored 
in the SDSS DR8 data base under the COLC and ROWC entries. The pixel 
coordinates in the SDSS DR8 data base correspond to the original detector 
coordinates because SDSS images are not corrected for optical distortions 
by resampling before o bject detection and photometry is performed on them 
I Stoughton et al. 2002). This implies that our calculated subpixel coordi- 
nates faithfully represent the subpixel coordinates of the detector pixels. 
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Figure 3. The r-waveband IPSMs as a function of distance from the centre 
of the map for each of the six detectors (red points). For clarity, the points in 
this plot corresponding to different detectors have been shifted along the x- 
axis by different fixed amounts. The master IPSM for all of the r waveband 
data is displayed as the image in the panel in the top left-hand corner with a 
linear colour-scale between =p0.25 mmag, and with the same orientation as 
the maps in Fig. fj] The corresponding radial profile of the master IPSM is 
plotted as the black points. 

the uncertainty in each magnitude offset. We note that the o val- 
ues are approximately the same (to within ~l-2%) for each sub- 
pixel region because the photometric observations are distributed 
uniformly over the pixel area. We force the mean of each image ar- 
ray to be zero, which is necessary for display purposes because of 
the general trend of the aperture magnitudes as a function of PSF 
FWHM and object brightness. The colour-scale in each of the im- 
ages is linear and runs between =f4<t and q=12a for the ugiz and r 
wavebands, respectively. 

The images above each plot in Fig. [2] may be interpreted as 
intrapixel photometric sensitivity maps (IPSMo for the regimes 
of under-sampled, normally-sampled, and well-sampled PSFs (left- 
hand, middle, and right-hand images, respectively) when measured 
by the SDSS pipeline aperture photometry routines with the 7.43" 
fixed aperture. These IPSMs show that for the ugiz wavebands we 
do not detect any systematic trends as a function of subpixel coor- 
dinates at the 4cr level (although there is a hint of a systematic trend 
for the i waveband), and they further highlight that the scatter in the 
magnitude offsets is greater than the corresponding uncertainties by 
a factor of ~1.7, which suggests the presence of some unmodelled 
systematic trends. We may translate the 4<r levels into upper limits 
on the systematic trends that may exist, and in Table [2] we report 
these 4<7 upper limits for the regimes of under-sampled, normally- 
sampled, and well-sampled PSFs for each detector in columns 5, 
6, and 7, respectively. In general, the upper limits that we de- 
rive are ~0. 1-0.3 mmag for the regimes of normally-sampled and 
well-sampled PSFs, and ~0.2-3.2 mmag for the regime of under- 
sampled PSFs. 

For the r waveband, the IPSMs in Fig.[2]reveal a smooth, ap- 
proximately radially-symmetric, systematic trend detected with a 

7 An IPSM jKavaldiiev & Ninko\Jll99Sh describes how the photometry of 
an unresolved object depends on the position of the PSF centroid within 
a pixel, and it is physically influenced by the pixel response function, the 
shape/width of the object PSF, and, if applicable, the TDI detector read- 
out mode. For under-sampled images, the effect is greatest, and as the PSF 
FWHM increases, the effect is washed out by the weaker PSF gradients 
over any one pixel. 
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Figure 4. Histograms of the magnitude offsets (mmag) for each filter that 
we calculate for use in calibrating the aperture magnitude measurements for 
the full set of selected observations from SDSS DR8 in Table [T] using data 
for the six camera columns combined. 



-~20<t range between the extreme values. The trend appears to be 
present for all six detectors in each PSF sampling regime. Hence, 
we have recalculated the IPSMs for each detector in the r waveband 
including all data regardless of the PSF FWHM. Then, in Fig. [3] for 
each detector, we plot the magnitude offsets from these IPSMs as a 
function of distance from the centre of the map (red points). Note 
that for clarity we have shifted the points corresponding to different 
detectors by different fixed amounts along the a>axis. Clearly the 
systematic trends are very similar between the six detectors, and 
therefore we have also derived a single master IPSM for all of the 
r waveband data, and we present this map in Fig. [3] as the image 
panel in the top left-hand corner. We also plot the radial profile of 
this master IPSM as the black points. Finally, we note that the form 
and amplitude of this IPSM is independent of right ascension, dec- 
lination, detector coordinates, object brightness, and the number of 
neighbouring objects within 15". 

The peak-to-peak amplitude of the r-waveband master IPSM 
is ~0.54 mmag, and it shows that objects observed with their cen- 
troid at the centre of a detector pixel are measured as ~0.54 mmag 
fainter than if their centroid is at the corner. This is in contrast 
to what we would expect if this systematic trend were due to the 
pixel response function (PRF), which generally results in objects 
being measured as brighter when their centroid is in the centre of 
a pixel compared to when t heir centroid is at the corner/edge (e.g. 
Kavaldjiev & Ninkov 1998). We also find that the detected system- 
atic trend as a function of subpixel coordinates is independent of 
the PSF FWHM, and therefore it cannot be the result of the PRF 
convolved with the object PSF. Furthermore, we do not observe 
the asymmetry we wo uld expect in the IPSM fo r detectors using 
the TDI readout mode dGibson&Hicksonlll992l) . This leads us to 
believe that the systematic trend is an effect introduced by the al- 
gorithms used to measure the aperture magnitudes, which is not 
so surprising when we find in Section |4~2l that the SDSS pipeline 
introduces a variety of systematic trends into the PSF magnitudes. 



4. 1.3 Correcting The Aperture Magnitudes 

Ideally, the systematic photometric trends that we have detected 
in the aperture magnitudes could be corrected by improving the 
image processing algorithms in the SDSS pipeline and then re- 
processing the imaging data. However, this solution is outside of 



the scope of this paper. Instead, we opt to correct the SDSS DR8 
aperture magnitudes by subtracting the fitted magnitude offsets Zk 
from the aperture magnitude measurements, which constitutes a 
post-processing self-calibration of the data. Since the correction of 
SDSS DR8 aperture magnitudes is of general interest to the as- 
tronomical community, we supply an IDL program (...weblink...) 
which can be used to correct a set of aperture magnitude measure- 
ments from SDSS DR8 for the systematic trends that we detected 
as a function of PSF FWHM and object brightness in Section |4.1.1| 
We do not include the corrections for the r-waveband master IPSM 
because of the very small amplitude of this effect. Note that these 
magnitude corrections should only be applied to aperture magni- 
tude measurements of point sources (since they were derived ex- 
clusively from such sources), and we stress that even though the 
subtraction of our fitted magnitude offsets from the aperture mag- 
nitudes constitutes a relative photometric calibration, any absolute 
calibration of the photometry may need to be redetermined. 

We may now assess the significance of our photometric cali- 
brations. We do this by first fitting the photometric model in Equa- 
tion Q] with no magnitude offset parameters to the aperture mag- 
nitude measurements in the full set of repeat observations in Ta- 
ble [2] and we calculate the associated chi-squared Xbcfore- We then 
subtract the magnitude offsets Zk corresponding to the detected 
systematic trends as a function of PSF FWHM and object bright- 
ness in Section 14.1.11 and which have an uncertainty of less than 
5 mmag, from the appropriate aperture magnitude measurements. 
Finally, we fit again the photometric model in Equation [T] with no 
magnitude offset parameters to the corrected aperture magnitude 
measurements and calculate the associated chi-squared Xafter- 

In Table [2] we tabulate the values of A\ 2 = Xbcfore — Xafter- 
By fitting and applying corrections to the aperture magnitude mea- 
surements, we have introduced N z ~ 90 new parameters between 
calculating the different \ 2 values, and we would therefore expect 
that A\ 2 ~ 90 if there are no systematic trends. The values of 
A\ 2 tabulated in Table[2]are much greater than 7V Z which indicates 
that our corrections to the repeat aperture magnitude measurements 
are highly significant. We note that the variance of \ 2 1S 2JVd f, 
where N&>t is the number of degrees of freedom. Hence, the signif- 
icance in units of a of the A\ 2 values is given by A\ 2 /V2N]m- 
For Xbeforo we have iVdof = AWa — Aobj (column 4 minus column 
3 of Table|2j, and for Xafter, we have iVdoi = A/data — N obi — N z « 
A^data - iVobj (sinc e N z < N obl ). In column 9 of Tabled we list the 
quantity A\ 2 /y2(JVdata — iVobj) for each detector. We conclude 
that our corrections to the repeat aperture magnitudes are signifi- 
cant at the ~20-40ct, 58-109<r, 74-120rr, 56-84cr, and 45-68<r level 
for the ugriz wavebands, respectively. We also note that if we re- 
produce Fig. [2] using the corrected aperture magnitude measure- 
ments, then the systematic trends as a function of PSF FWHM are 
not present any more, which confirms that applying the corrections 
has successfully calibrated the aperture magnitudes to compensate 
for this particular systematic trend. 

To conclude with our analysis of the SDSS aperture magni- 
tudes, in Fig.|4]we plot for each filter a histogram of the magnitude 
offsets (mmag) that we calculate for use in calibrating the aperture 
magnitude measurements for the full set of selected observations 
from SDSS DR8 in Tablefj] combining the data for the six camera 
columns. We see that there are a significant number of corrections 
that have absolute values >2 mmag and stretching up to ~ 10 mmag 
and beyond, where ~10 mmag is the point at which the distribution 
seems to flatten out in both directions (there are also some mag- 
nitude corrections that lie outside of the plot limits). These mag- 
nitude corrections are of the order of the quoted ~1% precision 
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Table 3. The properties of the selected repeat observations with PSF mag- 
nitude measurements from SDSS DR8 (see Section 13.3) for each of the 
30 SDSS imaging detectors. The data sets for each detector are organised 
by filter and camera column. In column 5, we report the amount by which 
the x 2 f° r the repeat observations decreases after calibrating the PSF mag- 
nitudes using the fitted magnitude offsets from Sections 14.2.11 & 14.2.21 In 
column 6, we list the Ax 2 values in units of a. 



4.2. 1 Trends As A Function Of PSF FWHM 

Fig-H]is the equivalent of Fig.[2]produced for the PSF magnitudes. 
The only difference in the plot schema is that in Fig. [5] the colour- 
scale in each of the IPSMs is in units of mmag rather than in units of 
a. For all wavebands, the colour-scale in these maps is linear, and 
it runs between q=3 mmag for the maps corresponding to the under- 
samp led PSF regime, and between =p0.8 mmag for the maps cone- 
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^gopding to the normally-sampled and well-sampled PSF regimes. 

Inspection of the plots of the fitted magnitude offsets as a func- 
tiun of PSF FWHM reveals that as the PSF FWHM increases from 



of the SDSS DR8 photometric measurements, and we believe that 
the application of these corrections to the aperture magnitude mea- 
surements will significantly improve the SDSS DR8 photometric 
precision in all filters. 



the critical sampling to ~2.2-2.5", the magnitude offsets "oscil- 
late" with an amplitude of ~2-6 mmag (the ug wavebands show 
the clearest examples) while also decreasing with a slight gradient 
(~l-5 mmag/arcsec; most notably in the g waveband). This be- 
haviour is largely independent of object brightness (except in the z 
waveband). It is also interesting to note that the form and amplitude 
of the trend is consistent between the different camera columns for 
each waveband. For PSF FWHMs greater than ~2.2-2.5", we find 
that as the PSF FWHM increases, the magnitude offsets increase 
sharply by up to ~ 20-50 mmag. For PSF FWHMs less than the 
critical sampling, the systematic trend tends to be that as the PSF 
FWHM decreases, the fainter an object is measured (except in the z 
waveband), with some evidence that this behaviour becomes more 
pronounced for even smaller PSF FWHMs (see the ri wavebands). 

We have tested that the detected systematic trends in the mag- 
nitude offsets as a function of PSF FWHM are independent of right 
ascension, declination, detector coordinates, subpixel coordinates, 
and the number of neighbouring objects within 15" by making 
similar plots to those displayed in Fig. [5] for various mutually- 
exclusive partitions of our photometric data. 

Clearly, the systematic trends in the magnitude offsets for 
the PSF magnitude measurements are more complicated than the 
trends observed in the magnitude offsets for the aperture magni- 
tude measurements, and this suggests that it is the algorithms used 
to measure the PSF magnitudes that introduce these trends into the 
photometry. The scientific literature contains only brief qualitative 
descripti ons of the PSF fitting procedures implement ed in the SDSS 
pipeline (jLupton et al.ll200ll ; lstoughton et al.ll2002h . and there are 
no published studies in which an analysis of the behaviour of these 
routines has been performed. Therefore, we do not attempt to of- 
fer an in-depth explanation for any of these systematic trends in 
terms of the algorithms used to measure the PSF magnitudes. We 
limit our scope to simply highlighting the existence and form of 
the systematic trends in the PSF magnitude measurements, and to 
speculating that the "oscillations" present in the magnitude offsets 
are probably due to some form of resonance between the PSF shape 
and the underlying pixel grid. 



4.2 PSF Magnitudes 

We now analyse the PSF magnitudes. We use the data correspond- 
ing to Table|3] Independently for each detector, we investigate how 
the PSF magnitudes correlate with various object/image properties. 
Using the same fitting procedure as that described in Section 14.11 
applied to the PSF magnitude measurements, we find that there are 
no significant trends in the fitted magnitude offsets as a function of 
detector row (Zk < 0.2 mmag) or column (Zk < 0.5 mmag), and, 
as with the aperture magnitude measurements, we do find a system- 
atic trend in the magnitude offsets as a function of PSF FWHM and 
object brightness. We also find a systematic trend in the magnitude 
offsets as a function of subpixel coordinates. 



4.2.2 Trends As A Function Of Subpixel Coordinates 

The IPSMs displayed in Fig. [5] show very clear systematic trends 
in each of the regimes of under-sampled, normally-sampled, and 
well-sampled PSFs. For the regime of under-sampled PSFs, we 
find that in the riz wavebands, objects with their centroid close 
to the centre of a detector pixel are measured by up to ~4-7 mmag 
brighter than if their centroid is at the corner of a detector pixel. We 
cannot discern this effect for the regime of under-sampled PSFs in 
the ug wavebands due to the lack of data in these wavebands with 
PSF FWHMs less than the critical sampling. Since these system- 
atic trends are not present in the regime of under-sampled PSFs 
for the aperture magnitude measurements, we can only conclude 
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Figure 5. Systematic trends in the fitted magnitude offsets Z k for the PSF magnitude measurements as a function of PSF FWHM, object brightness, and 
subpixel coordinates (see text in Sections |4.2. 1 1 & [ 4.2.21 for details). The data used to generate these plots are the selected repeat observations with PSF 
magnitude measurements from SDSS DR8 (see Section [33). The black, red, and green points correspond to PSF magnitude measurements brighter than 
16 mag, in the range 16-18 mag, and fainter than 18 mag, respectively. Asterisks represent magnitude offsets that fall outside of the plot range. The three 
images above each plot represent the IPSMs for the regimes of under- sampled, normally-sampled, and well-sampled PSFs (from left to right). The colour- 
scale is linear and runs between ^3 mmag for the under-sampled PSF regime, and between ^0.8 mmag for the normally-sampled and well-sampled PSF 
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Figure 6. The master IPSMs for each waveband as a function of distance 
from the centre of the map. These master IPSMs have been derived for each 
waveband using all of the PSF magnitude measurements over the six de- 
tectors that have corresponding PSF FWHMs greater than the critical sam- 
pling. For clarity, the points in this plot corresponding to different wave- 
bands have been plotted with different colours and shifted along the x-axis 
by different fixed amounts. The ^-waveband IPSM has a mean value of 
zero, and for comparison purposes, the maps for the other wavebands have 
had a constant added to the magnitude offsets so as to match the central 
magnitude offset between the maps. The images in the top panels display 
the master IPSMs for each waveband as a function of subpixel coordinates 
using a linear colour-scale between =p0.8 mmag, and with the same orien- 
tation as the maps in Fig. [5] 



that they are introduced into the PSF magnitude measurements by 
the algorithms used to measure the PSF magnitudes, and that the 
underlying cause is not a physical effect (i.e. a non-uniform PRF 
convolved with the object PSF). We speculate that a possible ex- 
planation for the trends could be a slight mismatch between the 
PSF model and the actual shape of the PSF. The fact that, for the 
regime of under-sampled PSFs, the exact form of the IPSMs varies 
between different wavebands and camera columns hints that the 
deficiencies in the PSF modelling depend on the exact PSF shape, 
which is unique to each detector. 

For PSF magnitude measurements with PSF FWHMs that 
are greater than the critical sampling, we find that for all wave- 
bands, the IPSMs in Fig. [5] have the same smooth and approxi- 
mately radially-symmetric form as the systematic trend detected 
for the aperture magnitudes in the r waveband in Section 14.1.21 
For each waveband, the trends are very similar in form and ampli- 
tude between the six detectors and the normally-sampled and well- 
sampled PSF regimes. Therefore we have derived a single master 
IPSM for each waveband using all of the PSF magnitude measure- 
ments over the six detectors that have corresponding PSF FWHMs 
greater than the critical sampling. We present these maps in Fig. [6] 
as the image panels along the top of the figure. In the same fig- 
ure, we also plot the radial profiles of the master IPSMs for each 
waveband. Finally, we note that the form and amplitude of these 
IPSMs are independent of right ascension, declination, detector co- 
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Figure 7. Histograms of the magnitude offsets (mmag) for each filter that 
we calculate for use in calibrating the PSF magnitude measurements for the 
full set of selected observations from SDSS DR8 in Table[T] using data for 
the six camera columns combined. 



ordinates, object brightness, and the number of neighbouring ob- 
jects within 15". 

The master IPSMs displayed in Fig. [6] are very similar be- 
tween wavebands with peak-to-peak amplitudes of ~ 1 .55, 1.61, 
1.68, 1.57, and 1.67 mmag for the ugriz wavebands, respectively. 
These amplitudes are approximately three times larger than the 
peak-to-peak amplitude of the same systematic trend detected for 
the r-waveband aperture magnitude measurements. Again, we be- 
lieve that this systematic trend is an effect introduced by the algo- 
rithms used to measure the PSF magnitudes. 



4.2.3 Correcting The PSF Magnitudes 



Our IDL program described in Section 14.1.31 (...weblink...) may 
also be used to correct a set of PSF magnitude measurements of 
point sources from SDSS DR8. The program subtracts the appro- 
priate fitted magnitude offsets Zu corresponding to the systematic 
trends that we detected as a function of PSF FWHM and object 
brightness in Section l4.2. 1 I from the PSF magnitude measurements 
to be corrected. Also, for each detector, we repeat the analysis from 
Section [47272] to derive two IPSMs, one for each of the regimes of 
PSF FWHMs less than and greater than the critical sampling. For 
each map, we force the mean of the magnitude offsets to be zero to 
minimise the effect on the absolute photometric calibration when 
they are applied. The program subtracts the appropriate magnitude 
offsets taken from these IPSMs from the PSF magnitude measure- 
ments to correct for the systematic trends that we detected as a 
function of subpixel coordinates in Section l4.2.2l For the ug wave- 
bands, we refrain from applying the magnitude offsets from the 
IPSMs calculated for the regime of under-sampled PSFs because 
these maps are very noisy due to the lack of appropriate data. 

In essence, the procedure for calibrating the PSF magnitudes 
is the same as for the aperture magnitudes except for the inclu- 
sion of the extra correction for the systematic trends as a function 
of subpixel coordinates, which introduces a further 50 parameters 
into the model for the photometric data. However, the number of 
extra parameters introduced for our corrections to the PSF magni- 
tude measurements is still insignificant compared to the total num- 
ber of free parameters. Therefore, we follow a similar analysis to 
that of the aperture magnitudes with regards to the significance of 
our corrections to the repeat PSF magnitude measurements, and in 
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columns 5 and 6 of Table [3] we tabulate the values of A\ 2 and 
A\ 2 1 \J '2(iVdata — iVobj) for each detector. We conclude that our 
corrections to the repeat PSF magnitude measurements are signifi- 
cant at the ~8-20ct, 35-132ct, 32-76<t, 28-68ct, and 9-2 1<t level for 
the ugriz wavebands, respectively. We also note that if we repro- 
duce Fig.[5]using the corrected PSF magnitude measurements, then 
the systematic trends as a function of PSF FWHM and subpixel co- 
ordinates are not present any more, which confirms that applying 
the corrections has successfully calibrated the PSF magnitudes to 
compensate for these particular systematic trends. 

In Fig. [7] we plot for each filter a histogram of the magnitude 
offsets (mmag) that we calculate for use in calibrating the PSF 
magnitude measurements for the full set of selected observations 
from SDSS DR8 in Table[T] combining the data for the six camera 
columns. Again, we see that there are a significant number of cor- 
rections that have absolute values >2 mmag and stretching up to 
~8 mmag. As with the aperture magnitudes, these corrections are 
of the order of the quoted ~1% precision of the SDSS DR8 photo- 
metric measurements, and we believe that the application of these 
corrections to the PSF magnitude measurements will significantly 
improve the SDSS DR8 photometric precision in all filters. 



5 COMPARING THE SDSS APERTURE AND PSF 
MAGNITUDE SCALES 

In this section, we investigate how the SDSS PSF magnitude scale 
relates to the SDSS aperture magnitude scale for the 7.43" fixed 
aperture. This investigation is important because the photometric 
calibration parameters derived by PAD08 are obtained by fitting 
the aperture magnitudes and the calibration is subsequently applied 
to the PSF magnitudes. This is a valid procedure for the calibration 
of the PSF magnitudes as long as they follow the same magnitude 
scale (to within the addition of an arbitrary constant) as the aper- 
ture magnitudes. Any deviation of the PSF magnitude scale from 
the aperture magnitude scale other than a simple offset renders the 
calibration parameters derived from the aperture magnitudes as in- 
compatible with the PSF magnitudes. Consequently, in such a case, 
the application of the calibration parameters to the PSF magnitudes 
will introduce undesirable systematic errors. 

We may use the selected observations from SDSS DR8 in Ta- 
ble Q] to compare the PSF magnitude measurements to the aper- 
ture magnitude measurements. We analyse the photometric data 
for each detector independently. Firstly, we drop the observations 
for which the SDSS pipeline failed to measure an aperture flux. 
In addition, we restrict our sample of observations to those ob- 
jects which have no neighbouring objects within 15". We do this 
to avoid the inclusion of objects with aperture magnitude measure- 
ments that have been contaminated by the flux of nearby objects, 
which would lead to artificially bright aperture magnitude measure- 
ments for these objects when compared to the corresponding PSF 
magnitude measurements. Next, for each observation, we calculate 
the magnitude difference Am = m ps f — m ap , where m ap and m ps f 
are the measured aperture and PSF magnitudes, respectively, cal- 
ibrated using the procedures described in Sections 14.1.31 & 14.2.31 
The uncertainty on Am is calculated as the square root of the sum 
in quadrature of the corresponding aperture and PSF magnitude 
measurement uncertainties. 

We partition the Am values into three groups based on their 
corresponding PSF FWHMs; namely, the under-sampled PSFs with 
PSF FWHMs less than the critical sampling of 0.93", the normally- 
sampled PSFs with PSF FWHMs greater than the critical sampling 



but less than 1.5", and the well-sampled PSFs with PSF FWHMs 
greater than 1.5". Then, for each grouping, we assign the Am val- 
ues to magnitude bins of width 0.2 mag using their corresponding 
PSF magnitudes, and we calculate the inverse-variance weighted 
mean of the Am values in each bin. In Fig. [8] for each detector, 
we plot these mean Am values as a function of PSF magnitude us- 
ing black, red, and green points for the regimes of under-sampled, 
normally-sampled, and well-sampled PSFs, respectively. In each 
plot, the error bars represent the uncertainties in the mean Am val- 
ues, and they are generally much smaller than the plot symbols 
(~0.03-0.2 mmag). We refrain from plotting mean Am values with 
uncertainties of greater than 5 mmag, and we mark mean Am val- 
ues that fall outside of the plot range with an asterisk of the relevant 
colour. The plot panels are arranged in the same grid layout as in 
Fig.H] 

If the PSF magnitude scale is the same as the aperture magni- 
tude scale except for a difference in photometric zero-points, then, 
independent of object brightness, the expected value of Am is the 
zero-point difference, and the plots of the mean Am values as a 
function of object brightness will be horizontal straight lines with 
a y-axis value equal to the zero-point difference. However, the fact 
that all of the plots in Fig. [8] show mean Am values that depend 
on the object brightness in a non-linear fashion (i.e. they do not 
form horizontal or inclined straight lines) indicates that the PSF 
magnitude scale is related to the aperture magnitude scale by a 
non-linear transformation. As a further complication, we may also 
see in Fig.[8]that this non-linear transformation is PSF FWHM de- 
pendent since the mean Am values show different behaviour for 
each of the regimes of under-sampled, normally-sampled, and well- 
sampled PSFs. 

We observe that, for any single waveband and PSF FWHM 
regime combination, the form and amplitude of the non-linear rela- 
tion between the aperture and PSF magnitude scales is very similar 
for all six detectors. For the griz wavebands, the fainter an ob- 
ject is, the brighter the PSF magnitude measurement is compared 
to the aperture magnitude measurement, where the difference be- 
tween the two types of magnitude measurements can reach up to 
~ 15-35 mmag for the faintest objects in our data set when observed 
under normal-seeing conditions. For the u waveband, we detect the 
opposite effect (except for camera column 2) with a difference of up 
to ~ 10-25 mmag between the aperture and PSF magnitudes for the 
faintest objects in our data set when observed under normal-seeing 
conditions. We also note that for the it waveband, the PSF magni- 
tudes are fainter than the aperture magnitudes by ~5-10 mmag, in- 
dicating that there is a zero-point difference between the two mag- 
nitude scales for this waveband. For the z waveband, the aperture 
and PSF magnitude scales diverge rapidly for objects with bright- 
nesses in the range 18-19mag. Also, for all wavebands, we find that 
the amplitude of the non-linearity in the relation between the aper- 
ture and PSF magnitude scales is greater for larger PSF FWHMs. 
Finally, we note that the erratic behaviour of the mean Am values 
at the bright-end for all wavebands is most likely a consequence of 
image saturation for the corresponding observations. 

We have tested that the detected non-linear relation between 
the aperture and PSF magnitude scales is independent of right 
ascension, declination, detector coordinates, and subpixel coordi- 
nates by making similar plots to those displayed in Fig. [8] for vari- 
ous mutually-exclusive partitions of our photometric data. Further- 
more, the plots in Fig. [8] are only marginally affected if we do not 
calibrate the aperture and PSF magnitudes using the procedures de- 
scribed in Sections |4. 1.31 & 14.2.31 This is to be expected since, by 
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SDSS imaging camera layout: Scan direction to the left 
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Figure 8. The mean difference between the PSF and aperture magnitude measurements as a function of object brightness (PSF magnitude). The data used 
to generate these plots are the selected observations from SDSS DR8 restricted to observations of objects that have no neighbouring objects within 15". The 
black, red, and green points correspond to the regimes of under-sampled, normally-sampled, and well-sampled PSFs, respectively. Asterisks represent mean 
differences that fall outside of the plot range. 



necessity, our magnitude corrections are calculated independently 
of object brightness. 

It is important to understand that our analysis may only deter- 
mine how the PSF magnitude scale relates to the aperture magni- 
tude scale. Due to the fact that there is a non-linear relation between 
the aperture and PSF magnitude scales, we may deduce that at least 
one of the scales must deviate from a pure asinh-magnitude scale. 
However, it is quite possible that both the aperture and PSF mag- 
nitude scales actually deviate from a pure asinh-magnitude scale, 
and that the non-linear relation that we find between the magnitude 



scales simply indicates that they differ from each other in a non- 
linear fashion. One way to investigate if a magnitude scale actually 
deviates from a pure asinh-magnitude scale is to simulate a set of 
images containing objects of known magnitudes, and then to pro- 
cess these images with the SDSS pipeline to extract the relevant 
magnitude measurements for comparison with the known magni- 
tudes. We note that both the aperture and PSF magnitude measure- 
ments that we have employed in this analysis have already been 
calibrated using the PAD08 modelling scheme. However, the cal- 
ibrations will not have affected the relation between the aperture 
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and PSF magnitude scales because, for each observation, the same 
calibration correction is applied to each type of magnitude mea- 
surement. 

We find that we can explain the main features of the non-linear 
relation between the aperture and PSF magnitude scales as, yet 
again, a consequence of sky over-subtraction in the SDSS images. 
The SDSS pipeline performs PSF photometry by scaling the image 
PSF to each object on the sky-subtracted images without fitting a 
local background (R. Lupton, private communication). Hence, the 
main difference between the SDSS aperture and PSF photometry is 
the pixel weighting scheme that is used when measuring the flux of 
an object. Aperture photometry weights all pixels equally within 
the aperture, whereas PSF photometry weights the pixels in the 
aperture using the image PSF, giving more weight to the image pix- 
els that contain most of the object flux. Hence, the effective number 
of pixels used in the PSF photometry measurement is less than the 
effective number of pixels used in the aperture photometry mea- 
surement, and so we would expect that if the sky background is 
over-subtracted from the SDSS images, then the PSF photometry 
measurement will be brighter than the aperture photometry mea- 
surement. Also, we have already seen in Section |4.1.1| that the aper- 
ture magnitude measurements of fainter objects suffer from larger 
systematic errors due to incorrect sky-subtraction, and so we would 
expect the PSF magnitudes to differ by a greater amount from the 
aperture magnitudes for fainter objects. These effects constitute the 
main features of the non-linear relation that we detect in Fig. [^be- 
tween the aperture and PSF magnitude scales. 

We note that the sky over-subtraction hypothesis cannot ex- 
plain the PSF FWHM dependence of the relation between the aper- 
ture and PSF magnitude scales. As the PSF FWHM increases, the 
PSF photometry pixel-weighting scheme tends towards the aper- 
ture photometry pixel-weighting scheme, and therefore, in the case 
of sky over-subtraction, we would expect the difference between 
the aperture and PSF magnitude measurements to decrease, which 
is the opposite of what we detect. 



6 SUMMARY AND DISCUSSION 

We have investigated the systematic trends in the SDSS DR8 pho- 
tometric data, extending the internal consistency checks performed 
by PAD08 on the aperture magnitude residuals from their photo- 
metric calibration model. From our analysis of the repeat observa- 
tions in SDSS DR8, we have discovered a systematic trend in the 
aperture magnitudes that is a function of PSF FWHM, the ampli- 
tude of which increases for fainter objects. This trend is present at 
the ~7-15 mmag and ~30-170 mmag level for the brightest (near 
the saturation limit) and faintest (19 mag) objects that we analysed, 
respectively, with the z waveband exhibiting the strongest trend. 
We have also detected a low-amplitude (~0.54 mmag) systematic 
trend in the r-waveband aperture magnitudes that is a function of 
subpixel coordinates. 

For the SDSS PSF magnitudes, we have discovered compli- 
cated systematic trends that are a function of PSF FWHM, object 
brightness, and subpixel coordinates. The trends that are a func- 
tion of PSF FWHM and object brightness are present at the ~10- 
20 mmag and ~ 20-50 mmag level for PSF FWHMs that are less 
than ~2.5" and greater than ~2.5", respectively. Again, it is the 
z-waveband PSF magnitudes that exhibit the strongest trends as 
a function of PSF FWHM and object brightness. The systematic 
trends in the PSF magnitudes that are a function of subpixel coordi- 
nates are clearly detected in all wavebands at the ~1.6 mmag level 



for PSF FWHMs greater than the critical sampling (0.93"), and in 
the riz wavebands at the ~4-7 mmag level for PSF FWHMs less 
than the critical sampling. In the ug wavebands, we cannot discern 
these trends as a function of subpixel coordinates for the regime of 
under-sampled PSFs due to the lack of appropriate data. 

To address the problem of these systematic trends, we have 
described a method for self-calibration of the SDSS photometric 
data which successfully removes the trends, and we provide an IDL 
program which can be used to calibrate a set of aperture and/or 
PSF magnitude measurements of point sources from SDSS DR8. 
However, it would be better to eliminate these trends by identifying 
the algorithm(s) in the SDSS pipeline that introduce them, and then 
improving the relevant algorithm(s) appropriately. Another way to 
avoid the need for a post-processing self-calibration is to include a 
set of appropriate terms in the PAD08 modelling scheme. 

We hypothesise that the sky over-subtraction problem dis- 
cussed by AIH1 1 is the main cause of the detected systematic trend 
in the aperture magnitudes as a function of PSF FWHM and ob- 
ject brightness, although there may also be some contribution to 
the trend from aperture flux losses. With regards to the complicated 
systematic trends that we have found in the PSF magnitudes as a 
function of PSF FWHM and object brightness, we do not speculate 
on exactly what is causing them because we do not have access to 
sufficient information on the details of the PSF photometry routines 
in the SDS S pipeline. Having s tudi ed the basic descriptions of these 
routines in Lupt on et al.l d200lh and Stought on et al. I d2002l) , we flag 
the following procedures as potentially containing algorithms that 
introduce the detected systematic trends into the PSF magnitudes: 

• The PSF model adopted for each image including the model 
for the spatially variable terms. 

• The procedure for fitting the PSF model to each object in or- 
der to measure the PSF flux, which includes sinc-resampling of the 
image pixel data for the objects. 

• The procedure for determining the aperture correction that is 
applied to the PSF magnitude measurements. 

The systematic trends in the PSF magnitudes as a function 
of subpixel coordinates for the regime of under-sampled PSFs are 
also most likely caused by some aspect(s) of the above procedures, 
and, specifically, we suspect that they are due to a slight mismatch 
between the PSF model and the actual shape of the image PSF. We 
are unable to offer an explanation for the unexpected systematic 
trend as a function of subpixel coordinates for normally-sampled 
and well-sampled PSFs that is present at the same amplitude in all 
wavebands for the PSF magnitudes, and which is detected in the r 
waveband for the aperture magnitudes. 

We also found that the SDSS aperture and PSF magnitude 
scales are related by a non-linear transformation that departs from 
linearity by ~l-4%, implying that at least one of the scales departs 
from a pure asinh-magnitude scale. Furthermore, this discovery in- 
validates to some extent the application of the PAD08 calibrations 
to the PSF magnitudes, since the calibration model parameters are 
derived from fitting the aperture magnitudes. To avoid the problems 
introduced by the non-linear relation between the magnitude scales, 
it would be necessary to either transform the aperture magnitudes to 
the PSF magnitude scale before performing the calibration model 
fit, or to transform the PSF magnitudes to the aperture magnitude 
scale before applying the photometric calibrations derived from the 
aperture magnitudes, or, better still, to simply identify and rem- 
edy the cause of the non-linear relation in the first place. We have 
argued that the main cause of the non-linear relation between the 
magnitude scales is most likely yet another consequence of sky 
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over-subtraction in the SDSS images. Finally, we report the detec- 
tion of a ~0.5-l% zero-point difference between the aperture and 
PSF magnitude scales for the u waveband. 

Our results indicate that there is still room for improvement 
in the relative photometric calibration of the SDSS photometric 
data. PAD08 report that they detect systematic trends in the resid- 
uals of the SDSS aperture magnitudes for the PAD08 model "at 
the ~0.5% level". They also perform simulations of the data that 
include random walks in time of the atmospheric extinction coef- 
ficients to simulate temporal atmospheric variations, and they fit 
the simulated data with the photometric calibration model to see if 
they can reproduce the systematic errors at the ~0.5% level that 
they have detected in the real data. PAD08 find that their simula- 
tions produce systematic trends at the level of ~13, 8, 8, 7, and 
8 mmag in the ugriz wavebands, respectively, and therefore, they 
hypothesise that it is the unmodelled temporal atmospheric varia- 
tions in the PAD08 modelling scheme that are the cause of the re- 
maining systematic trends in the SDSS photometry. However, the 
results from their simulations depend heavily on the assumptions 
about the temporal atmospheric variations. We have shown, in con- 
trast, that there are still systematic trends not related to temporal 
atmospheric variations in the SDSS photometry at the 0.5% level 
and above for both the aperture and PSF magnitudes. Hence, it is 
possible to further improve the precision of the relative calibration 
of the SDSS photometric data, which, if implemented, will be of 
great benefit to the astronomical community. 

The SDSS photometry is now routinely used to calibrate 
other photometric data, inc luding other large scale surveys (e.g. 
Brown et al. 201 1, Sesar et al Will etc.). Specifically for this pur- 
pose. Ilvezic et alj d2007l) have produced a standard star catalogue 
for Stripe 82 based on the SDSS repeat observations. Therefore, 
the correction of the systematic trends and magnitude scale non- 
linearities that we have found in order to produce the most precise 
SDSS photometric data is crucial. 

Finally, all of the issues that we have highlighted with the 
SDSS photometry are likely be releva nt to ot her upcoming s urveys 
( e.g. PanSTARRS -iKaiser et alj200l DES - lFlaughej20ol, LSST 
- Ilvezic etai]|2008l etc.) in their analysis of the photometric preci- 
sion that they achieve. Our paper emphasises the fact that the rela- 
tive calibration of the well-established SDSS photometry could still 
be improved further for little effort, and that a thorough investiga- 
tion and understanding of the systematic trends that are present in 
the photometry is very important when constructing a photometric 
calibration model so as to include all of the relevant terms. 
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